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CORRECTIONS TO THE THIRD VIRIAL COEFFICIENT FOR THE 
MODIFIED BUCKINGHAM EXPONENTIAL-6 POTENTIAL FOR HELIUM 


by 


John C. Edwards* 


ABSTRACT 

The objective of this work is to determine the contribution of 
nonadditive interatomic interactions (attractive and repulsive) to 
the classical third virial coefficient for an assumed modified 
Buckingham Exponential-6 pair potential for helium. Methods of 
statistical physics were used to determine corrections to the third 
virial coefficient. In these calculations the results of Kihara 
(4)* for the nonadditive three body attractive interaction between 
three spherically symmetric atoms, and Sherwood, De Rocco, and 
Mason's (10) asymptotic Dadi pees to Jansen's (3) expression 
for nonadditive three body repulsive interactions based upon an 
assumed Gaussian charge distribution were employed. Computer 
calculations were executed for several reduced temperatures. 
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Underlined numbers in parentheses refer to items in the list of 
references at the end of this report. 
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These calculations were carried out using programs HA 2983, HA 2860, 
and HA 3068, which were prepared by the Automatic Data Processing 
Branch of the Helium Activity. The results, based on improved 
numerical calculations, arsenide the conclusion of Sherwood and 
coworkers (10) and Sherwood and Prausnitz (11) that the nonadditive 
effects tend to cancel each other, and hence the total correction 


is not strongly dependent upon the pair potential selected. 


INTRODUCTION 

As part of Physics Project 7017 at the Helium Research Center, 
consideration is being given to a determination of the correction to 
the third virial coefficient due to nonadditive interatomic inter- 
actions for helium. The work of other investigators was followed. 
The total interaction among three atoms was written as the sum of 
the pair interactions, and nonadditive corrections to the attractive 
and repulsive parts of the pair interaction. The pair interaction 
selected was the Modified Buckingham Exponential-6 potential 
(hereafter referred to as MBE6). This selection was based on 
its functional resemblance to the theoretical potentials for helium 
of Slater and Kirkwood (12) and Margenau (5). The short range 
repulsive part of a typical pair potential is based upon first 
order perturbation theory with exchange, whereas the long range 
attractive part of the pair interaction is based upon second order 
perturbation theory using wave functions that do not satisfy the 
Pauli principle. The nonadditive repulsive correction considered 


here comes from first order perturbation theory with exchange applied 
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to three atoms, where each atom is assumed to be surrounded by an 
"effective" electron represented by a Gaussian charge distribution, 
the model of Jansen (3). The nonadditive attractive interaction 
considered comes from third order perturbation theory applied to 
the dipole-dipole interaction among three atoms, that neglects 
overlap and exchange. This model was developed by Axilrod and 
Teller (1) and Kihara (4). These results, and Sherwood'’s approxi- 
mation to Jansen's result, were used to calculate the linear 
corrections to the third virial coefficient, as well as the 
additive part of the third virial coefficient, for several reduced - 
temperatures. 

Corrections to the third virial for the Lennard-Jones (12, 6) 
potential have been calculated by Sherwood and coworkers (10). 
Sherwood and Prausnitz (11) have computed the correction to the third 
virial due to nonadditive attractive interactions for the MBE6 


potential for vy eh ip & This work extends their results by considering 


For definition of y see representation of MBE6 potential on page 10. 
both attractive and repulsive corrections to the third virial for an 
assumed two body MBE6 potential that has y set equal to 12.4 on the 
basis of Mason and Rice's (6) analysis of helium data. 

Although the cumulative correction to the third virial coefficient 
is small, it is important to have as much information as possible for 
helium on the relationship between the interatomic potential and the 


third virial coefficient in order to obtain the best equation of 
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state, and a description of the relationship of the microscopic 


interactions to the observed macroscopic thermodynamic properties. 
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NONADDITIVE ATTRACTIVE INTERACTION 
Second order interactions between spherical neutral atoms are 
additive if the Pauli principle is not taken into account in the 
wave functions (5, 8). The result is the usual additive, -~ (atomic 
ee ae term. The next consideration is third order pertur- 
bation theory applied to dipole-dipole interactions among three 
atoms, and yields the nonadditive result in the absence of overlap 


and exchange 


eee : 
Ag = Grater ) [2 +3 cos 6, cos @, cos 6, | (4) (1) 


where 8. is the angle formed by the interatomic distances 4 and 
“Fee 


Se, 


ik’ ’ The quantity y is related to the static polarizability and 
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the constant » in the second order attractive term -4 r, where r 


is the interatomic separation, by the relationship 
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to express A Patt as 


Perey 
A =z - 3, He _att \ i (x : te 
Vatteiders wa cb 12° 13°23 
ij 
x [1 To3 COS o, cos 8, cos @|- (11) 


The quantity A Pore is positive if all o. are less than 117°, and 
consequently weakens the attractive forces. If one of the 0. is 
greater than 126°, then A : Pipa is negative. However, the latter 
case is less probable. 

It should be noted that Graben, Present, and McCulloch (2) 
considered the second-order single-overlap exchange contribution. 
However, aS pointed out by Sherwood and coworkers (10), it is 
inconsistent to consider such a correction without including the 


corresponding term in the pair potential. 


NONADDITIVE REPULSIVE INTERACT LOD 


As mentioned earlier, the repulsive part of the pair potential 


(3) 


is based upon first order perturbation theory with exchange. Extending 


the analysis to three helium atoms, Rosen (9), used hydrogenic wave 


functions and considered interatomic and intra-atomic exchange of 


electrons in a first order calculation, to determine the nonadditive 
repulsive interaction. The results are in the form of a number of 
cumbersome quadratures. In order to have a simple analytic form for 


Free Sherwood and coworkers (10) approximated the Gaussian model 
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of Jansen (3) for near-equilateral triangle configuration of atoms 
with fairly large interatomic separation. It is in this spirit that 
the calculations of this report were pursued. 

Jansen (3) bpp epee each atom to be surrounded by a Gaussian 


charge distribution, namely, 
p(x) = (BM/m)> exp ( -87r”) (4) 


which is equivalent to an "effective" electron on that atom, The 
distance from the center of the given atom is denoted by r, and 8 
is a constant dependent upon the species under consideration. 
Single interatomic exchange is accounted for between two and three 
atoms. Jansen used the wave function associated with the above 
charge distribution, and constructed the properly antisymmetrized 
wave function for three atoms, to calculate the nonadditive cor- 
rection to the repulsive interaction. Sherwood's approximation 


(see above) yields, for three atoms, 


Jj = 
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and for a pair of atoms the approximation reads 
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A.. = exp (| -B r,./4); r,. is the distance between atoms i and j, 
1j 1j ij 
Yi (jk) is the distance from atom i to the midpoint of the line joining 


atoms j and k, and e is the electronic charge expressed in esu. 
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Equations 5-a and 5-b imply 


1/2 1/2 
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Equations 3 and 6 are expressed in terms of an assumed two 
body potential. The explicit form of the MBE6 potential used in 


this report is 


a { oo Ce 
p(r) = 4 E 
ro sl$ oe [r- 2) -(2"] «> 


where ¢ is the depth of the potential well at its minimum rym 
locates the position of the spurious maximum, and y influences the 
steepness of the repulsive part of the potential. 

In table 1 a listing of a tcep and A? ott for an equilateral 
triangle configuration and the MBE6 pair potential » is given for 
1A<r<5 A. In the MBE6 potential the force constants determined 
by Mason and Rice (6) from an analysis of second virial coefficients 


and viscosity data for helium were used. These are 


tam 12 .4 

e/k = 9.16° K 
r = 3.135 A’. 
mh 


o3 
For the atomic polarizability the value ~ = .201 4 reported by 


Kihara (4) for helium was used. 
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TABLE 1. - Computer results for p(r), A Pee? and Ag. for 
ee ee eee ow eee Cp 


to, cm 


1.CO00000000E-08 
1-199900Q000E-08 
1. 2000000000E-08 
© 20)0000000E-08 
-“CO0000000E-08 
© 5O000000000E-08 
1 -GOONND00000E-08 
1-7000Q000000E-08 
1 -FONND0D00NE-08 
1 -°N000N0000E-08 
2 -CNDAD00000E-08 
2 e- LNDNADHNVOVE-08 
2 .2ONNNN0000E-08 
2 eB NNDDNQ0000E-08 
2 .4000000000E-08 
2 SO0NNN0000E-08 
2 -400NN00000E-08 
2 6 7ONN000000E-08 
2 -P000000000E-08 
2 -S00000000NE-08 
3 .OCONDN0LO0E-98 
3..1L6000N0000E-08 
“3 .2000000000E-08 
3 .300NN00000E-08 
3.4NNNNND000E-08 
3 5000000000E-08 
3.6000000000E-08 
3 7NONDN00N00E-08 
3. AOCONDOD00E-08 
3. 9D0MDD0000E=-08 
4&.0009000000E-08 
4.1000000000E-08 
4.2900000000E-08 
4.3000000000E-08 
4.4000000000E-08 
4 S50000N0000E-08 
46.6000000000E-08 
4.7909000000E-08 


et te 


4. 8000000000E-08 © 


4.9000000000E-08 
5 -Q000000000E-08 


te Nay a acy ices se ng 


- Patt? <*8* 


48213811910E-13 
2 .0447362804E-13 
923441597551 E-14 
4.5465429545E=14 
2 3335609475E-14 


“1.2541518924E-14 


7-0160330376E-15 
4.0656611227E=-15 
2 64306303889E=-15 
1.4941 323214E-15 
9.4167601388E-16 
6 .0701275473E=-16 
329936255476E-16 
2.6768309506E=-16 
1.8250312021E-16 
1.2638961509E-16 
8.8799667080E-17 
6 .3226274407E-17 
4.5577362256E-17 
3.3234525930F-17 
2 24495154148E-17 
1.8235439076F-17 
1.3703189526E-17 
1.0388 336536E-17 
729407443804E-18 
621172900103E-18 
407473249785E-18 
3.7098466046E-18 
269182271910E-18 
2-3098831248E-18 
1.8392109646E-18 
-1.4727083826E-18 
1.1855717865E-18 
9.5930204236E-19 
7.8000498977E=19 
663717517269E=-19 
5.62281854505E-19 
4.3081493239E=-19 
3 25645140667E=-19 
22960784) 627E-19 
26468 5471698EW19 


he aed ergs 


-2250154248612E-12 
my 4eIS01 632982 
-823649514852E-13 
-4.8103529178E=13 
~2-01980479053E-13 
-1.5773054884E-13 
-“9.0004171142E-14 


=5,12Z57153315E=14 — 


—2-9140957T600E-14 
-1.6541573821E-14 


—923766452454E-15 .- 
eA US DISS Gag Ee = 15.5 


-3.0019865966E=-15 
-1.6958748006E=-15 
-9.5712343570E-16 
-5.3971538020E=-16 
-3.0409825246E~-16 
-1.7121491438E-16 
-9.6332135330E-17 
-5.4165609059E=-17 
-3.0438112635E-17 
-1.7095050963E=-17 
-9.5961494402E-18 
—5.3840792510E-18 
-3.0194399227E-18 
—1.6925964166E-18 
9 64842524582E-19 
-5.3123316212E-19 
-2.9744630254E=19 
-1.6648746492E=-19 
—9 6 3156184654E-20 
-5.2108201787E=-20 
-2.9138773229E-20 
-1.6289708271E-20 
-9.1041179914E-21 
—5 (186865523 7E-21 
26641550483 7E-21 


-1.5869302906E=21 , 


—8e8605769194E-22 
-4.94620Z23 7719E=22 
22 7T605225844E=22 


equilateral triangle configuration for - 


o(r), ergs 


3.1869364899E-12 
263990061237E=12 
1.7203737635E-12 
1.2009687471E-12 
8.2418732067E=13 
5 056738 86962E-13 
3.7506377186E=13 


284952 50S099E~1 3% 


12645057688 3E=-13 
120737041265E=13 
6.9240594459E=-14 
4229717949335 =-14 
207352611761E-14 
126517882923E-14 
925298497 378E=15 
5 «0846209 104E=-15 
2 6 3089338536E-15 
6 62064057220E-16 
-3 .6644972013E-16 
-9.0714697787E-16 
-1.1685202121E=-15 
—1.2592803542E=-15 
-1.2494962167E-15 
-1.1837795386E-15 
-1.0900718910E-15 
-9.8547856490E-16 
-8.8012552442E-16 
-7 2796968539 1E-16 
-6 «8709432924E-16 
-6.0351502025t-16 
—5.2914466758E=-16 
—-4.6359860954E=-16 
=4.0619780389E-16 
—3 .5613790255E-16 
-3.21258959448E-16 
—2.7475529322E-16 
—2.4189852865E-16 
—2.1335663808E=-16 
-1.8854355650E=16 
1.669470 1 807E=16 
-1.4812286043E=-16 
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ft is seen from table 1 that the individual corrections to the 
three body interaction fall off rapidly with increasing interatomic 
separation, and that A Hep dominates over A eee up to r= 3 ie at 
which point A 7. becomes the dominant term. The result is that, 
except for a small region past the zero of the potential, the total 
correction has a quenching influence on the two body potential. 

For the equilateral triangle case, Rosen's (9) quantum 
mechanical calculation for helium yields as an asymptotic approxi-~ 


mation for r/r ral is 


= Prep rs xr 
sar -1.15 exp -5 -86516 e (8) 
rep 


This quantity is compared in figure 1 with ——S2 for Sherwood's 
Prep 


Figure 1. - Ratio of First Order Nonadditive Repulsion Energy to 
Total Pair Repulsion Energy for an Equilateral Triangle 


Configuration. 


approximation to Jansen's result for equilateral triangle configuration, 


namely 


tans 2(4-1) ow[-(Ge 2) ° 


y Prep 


and with equation 6 without the Gaussian dependence, that is 
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FIGURE |.-— Ratio of First Order Nonadditive Repulsion Energy to Total 
Pair Repulsion Energy for an Equilateral Triangle Configuration. 
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Freeh. | 223 @ - /3) \retew Ev - 3) « E (10) 
} a m m 


Prep 


A least squares procedure fitting equation 9 to equation 10 
over the range 0.5 arma 1.8 yielded § = 4.935575. It is seen 
from figure 1 that faa is very good agreement between Rosen's 
quantum mechanical calculation and Jansen's model modified with the 
assumed pair potential. The explicit Gaussian-dependent model is 
in poor agreement with either of these cases. This is expected 
since equations 8 and 10 are of similar functional form except for 


the factor /r in equation 10, whereas equation 9 involves a quadratic 


argument in the exponent. 


EXPRESSION FOR THE THIRD VIRIAL COEFFICIENT 
In the treatment of Mayer and Mayer (7) the classical third 


virial coefficient can be expressed as 
2 " 
C= (4b, B 2b.) N (11) 


where by is the k-th cluster integral. The cluster integrals are 
related to the microscopic regime through the interatomic potential 


as follows 


bo = a \ | exp (-8 $ (4p) 1 dr,dr, 


= 2n | | exp (-8 @. (ry) - 1 pda (12-a) 
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(1 . a — =) vie a 
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— 2 &* 


(e-8h) 


mS SS [owe 4 webs) «om (24 ch) 


i | + 9 + 4 >+,7,° 
- exp (-6 $ (x3 .r3)) - exp (-8 $ (<t.3)) eM drjdrjdr, 


' i ( \ \ | exp (- ° pet 3eF a1), uty (+2 "19)) 
- exp (-8 ) (r54)) - exp (-8 $ (r43)) 


+ 

2| ¥ 12723731 97 12°F 299731 . 

where 8 = 1/kT, k is Boltzmann's constant, and T is the temperature. 
Decomposing (249 2b 32631) into its additive and its nonadditive 


contributions, 
3- 
= + + 
(049 28 p5 91) if PSEA as tae Oe 
{£5 = 1 


and introducing the quantity 


fi = exp (- Cr) (x, ,)/kz) - 1 


it is possible to express C(T) as 


C(I) = C44 () +ac(z) (10) 


where 


5 | 
Bay? 6 
Caaa‘?) 3 \ \ Fioto3f 13% 1271372397 1297939 13 


14 


(12-b) 


(13) 


(14) 


(15) 


(16) 


& oy 7). % te) m9 “ ees A 
cress fe ee : eon «(Gate ® Jo ’ 
ea i 

Go ?. *) qxe = oateapee a ll De: 


“Ge. 9. a) aca Mes : ee + 


w 
ii, 
va 
a =) 


is ae 


fast To wer : 


) E ‘ 
* 7 | : 
(£4) aan? & on? a+ Ga) & BY, - Ceasegt2r 9 oa’ 


; Raat (2908 ; 
(at) , | (mA\4 2) y ) a te . 
cat) @ Mat & om qe 


(ar) artery aby worst wuts 7 a. 


4 é 5 # in te a7 


15 


AC(T) = a \ \exp e Le &p weet) 


— a L | ae “dra 
| exp ( ur) \ ) "1913230" 120° 13, °23° (17) 
If it is assumed further that Ag is a small quantity, then 


29 ” 
a(t) as SSE p \ \ \3 bel (-2 Le )) Fyoh13%Q 99 yo IT 14893 (18) 
i<j 


In this linear approximation, the quantity AC(T) is decoupled as 


= + 
Ac = AC + ACL (19) 


where AC is a function of the assumed pair potential, the 


rep (att) 
temperature, and A ree (a Part) ° 


Sherwood (10, 11) in effect rewrote equations 17-19 in terms 


of reduced units as 


£(TR) = Ce (TROT ACK (TH) +: ACK. (TX | 20 
CH(TH) = Chi (TH) + AGH (TH) + ACH (re) (20) 
C : : , ee 2 
where C* = “> and Re is the covolume, which is equal to 3 7 Ni 
ir 
kT o4 
He introduced the reduced units r* = oat ce yas. p* = z a* a 
m o 
ex = ee and changed the integration pattern by the transformation 
e 
r¥, = xr*, 73 amie a 33 = r* to obtain 


£ (e*) £(yr*) £(xr*) xyr®’dxdydr* = (21) 
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where f(s) = exp (- o*(s)/T*) - 1, 


3 
* 4 foe) 
ates eo) (2) 
a* T* \o 


be4 ptartay oa ers 
2 *4 
eee 
(xy) “r* [1 + 3 cos 8, cos 8, cos 2, | dxdydr* baz) 
13 
. 2 x Le 
Ac* r Al Ay o* (r¥.) 
LEED ag Ue i . 
jew ene ee) T* _ \ \ mis © T* )( * Crp de* GF 4)9* (x§)) 
= y rep rep rep 
a ert 


2 2 -1 -1 =] -1 =] 
ke ~e + a ae Si 
(xy) x lea SA, do (13) d3 (12) a y 1| gaa ae (23) 
where 


cae 
Asya etry 


d =< \I 2 i 
21135) fa) Va he che 


ae 
43012) “72 \¥? +4 -ig 


The point x = > poses mathematical difficulties in computing 
soy. Sherwood and coworkers (10) avoid this region by noting 
that the asymptotic expansion of the Gaussian model, equation (5-a), 
fails for any angle greater than 100°, and that the full Gaussian 


model is near zero in this region. Hence, the integrand was set 


equal to zero for qd) (23) < 0.4196. 
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The temperature range for the above expressions is valid for 
T >> 75° K, since for lower temperatures the thermal wavelength for 
helium is greater than the molecular diameter, and consequently 
quantum effects become important. 

The quantities Cr ad 
IBM 1620 computer* for T* = 10, 20, 30, 40, and 100 using Simpson's 


» AC* //e*, AC* /a* were computed on an 
rep att 


* Reference to specific models of equipment is made for identi- 


fication only and does not imply endorsement by the Bureau of 


Mines. 


rule for numerical integration. The programs prepared by ADP were 


HA 2983, HA 2860, and HA 3068, respectively. Table 2 gives the 
results of these calculations for ACE Ve» corresponding to 
three cutoff angles considered: 90°, 100°, and 110°. The only 
parameter specified was y, and Mason and Rice's (6) value of 12.4 


was used. 


TABLE 2. = Computer results for Cr aa? O en Ve*> and ACs _/a* 


for five reduced Saar 


AG* Skis 
ACh Vvé 


T* 


| 10 | -0.27 ae 32288 31525 

| 20[ =,2228607 |= 25651750 | ; 

sof aes | —ceecseia| cateonree | asses cuesaaes 
= 2195 S047 = ,22921635 | |e ai 
100 | -.17734976 | -.21233969 
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COMPUTATIONS OF C*¥ 
add 


The integration scheme for Cr ad was Ar* = .0125 for 0 < r* < 3 
and Ar*® = .025 for 3 < r* < 6 and AY = AX = .0125 for 0 < r* < 6. 
This is a refinement pn the interval size selected by Sherwood 

and Prausnitz (11). The error involved in using Simpson's rule 

is estimated to be within 1 percent by comparison of results for 
successive interval sizes. The discontinuities which occur in 

the integrand whenever r*, xr*, or yr* equals m* were accounted for 
by applying Simpson's rule only to continuous regions, avoiding 
integration over discontinuities. The computer program would 

adjust for the discontinuities by selecting the next smallest sub- 
interval size for the integration according to the first larger 

even number of subintervals which would produce a subinterval smaller 
than (or equal to) the specified subinterval size. It seems to this 
author that this was not considered by Sherwood and Prausnitg (11) 
for the following reason. Their value of 0.2798 for Ce ag forsy «IS 
and T* = 10 was reproduced using their subdivision scheme. A result 
of 0.2837 was obtained by using the same subinterval size, but 
accounting for the discontinuities. This seems to be a significant 
computational improvement. This aspect was incorporated into the 


computations of mone and AC + also. 


COMPUTATIONS OF ACE OM e* 
The integration pattern selected for Ack. We* was Ar*® = .02 
for 0 < r* < 2.5986383 and AY = AX = .015625. As mentioned above, 


discontinuities were accounted for in the integration. The results 
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are larger in absolute value than those found by Sherwood and co- 
workers (10) for the IJ (12, 6) potential. It is seen from table 2 
that a change in eutoff angle from 100° to 90° or 110° produces 
changes in Ace Me* eon 7.5 percent to 16.5 percent, where the 
discrepancy increases with increasing temperature. However, the 
Gaussian model itself is subject to criticism, as noted by Sherwood 
and coworkers (10). The net effect of the contribution of Ack We* 
is to lower the third virial coefficient. The same procedure was 
used to select the best interval size as above, which is within an 


error of 0.4 percent. 


COMPUTATIONS OF AC*_, /a* 
att 
The integration pattern for Ace /a* was Or« = .0l for 0< r*< 2, 
Ar® = .02 for 2 < r* < 10, and AY * AX = .015625. The subinterval 
length was one-half that used by Sherwood and Prausnitz (11). The 


effect of discontinuities were accounted for as before. The net 
effect of this correction to the third virial is to raise the curve. 
The error in the subinterval size used involves no error greater 
than five one-thousands of 1 percent. 


In all cases the ecutcff error in the r* integration was 


estimated to be within oneetenth of 1 percent. 


Ts > Ck * d AC* 
COMPARISON OF a a’ pecape an Cr tt 


d 


Table 3 lists hand calculator results applied to table 2 for 
specified values of ¢ and aw for heliwn, and the quantity aor eP is 
listed for Mason and Rice's (6) value of é. Ack is listed for 


two extreme cases of a: one is the value in the literature by 
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Kihara (4), and the other is the theoretical value calculated by 
Slater and Kirkwood (12). It is felt that the actual value should 
lie somewhere in between the two extremes. The results of table 3 
show that the contribution AC increases with a. 

The a cia as Lo Shr and ACS et are comparable in magnitude but 
of opposite sign, which, as reported earlier (10, 11) yields a small 
net contribution. In table 4 it is seen that |ac| is less than 


2 percent of eats 


s) 


ae ae 
oe @x3 = = * 
TABLE 3. ACT ee for a 201 A anda eh ily ACK 


for ¢/K = 9.16° K, and AC* 


AC® AC* 


“Kihara (4). 
= Slater and Kirkwood (12). 


TABLE 4. = Percent ratio of JAc| to Ics aal 
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DISTRIBUTION FUNCTIONS FOR AC*, AC* , and AC* 
’ rep att 


The quantities AC* and AC*,_. can be written as 
rep att 


Cr ep = Jex n\ \ \ f£(x,y,r*) dxdydr* (24-a) 


hear arr Ht \ \ e(x,y,c*) dxdydr* (24-b) 


where N and M are constants 
and 


ACe = AC | + AC* (24-c) 
act rep 


The normalized distribution functions which have been averaged 
over all configurations for some given separation r* of two of the 
three atoms are 

N £(x,y,r*) dxdy 


Ce (r%) = 7p yee 
rep. || | \ \ \ f(x,y,r*) dxdydr*| 


M g(x,y.t*) dxdy eons 
Co. Az*) & f— 25=-b 
ae || | \ \ \ g(x,y,r*) dxdydr*| 


sat grt 
rep ¥ o a lo 
+ * * 25- 
Jac*| Coan [Ac*] a ) rae 


C* (r*) 


The values of C* , C¥ 
re 


ae aces and C* as given by equations 25-a, 25-b, 


23 
and 25-c are plotted in figures 2 and 3 for @ = .201 A” and Mason 


7, 


Figure 2. - Distribution Functions for AC*, 3 Sao and AC* at T* 


100. 


Figure 3. - Distribution Functions for AC*, OS oe and ACT et at T* 
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FIGURE 2.-— Distribution Functions for Ac”, AY eB 
and Corts at T= lO. 
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and Rice's (6) force constants for T* = 10 and 100 respectively, 

with a cutoff angle of 100°. These quantities were readily avail- 
able from the computer output since the integrand of the r* 
integration, the result of the integration over x and y, ie printed 
out for each value of r* considered. It is seen from the graphs that 
C*(r*) has more positive (negative) contribution for T* = 10 (100), 


which corresponds to the signs of AC* for these two cases. 


CONCLUS IONS 

The expected result of a near cancellation occurred for the two 
corrections to the third virial coefficient. It was felt necessary 
to explore the full effect of these corrections for an assumed MBE6 
potential because it is more realistic than a Lennard-Jones (12,6) 
potential for heliun. 

It seems that greater precision in computation has been achieved. 
However, it still remains an unanswered question as to why effects 
arising from two independent physical phenomena should nearly cancel. 

The computations presented here (T* > 10) correspond to the 
high temperature tail of the third virial coefficient for helium. In 
the very low temperature range the deviation, AC, becomes quite 
significant. Extending the calculations down to this range would therefore 
involve a quantum mechanical calculation of the second and third cluster 


integrals. I hope to return to these considerations in future work. 
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